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I analytically study the plateau of the magnetization curve at M/Ms = 1/3 (where Ms is the saturation 
magnetization) of the one-dimensional 5 = 1/2 trimerized Heisenberg spin system with ferromagnetic (Jf)- 
ferromagnetic ( Jp)-antiferromagnetic (J a) interactions at T = 0. I use the bosonization technique for the 
fermion representation of the spin Hamiltonian through the Jordan- Wigner transformation. The plateau 
appears when 7 = Jf/Ja < 7c, an d vanishes when 7 > 7c, where the critical value 7c is estimated as 
7c = 5 ~ 6. The behavior of the width of the plateau near 7c is of the Kosterlitz-Thouless type. The 
present theory well explains the numerical result by Hida. 
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Recently Hida 1 numerically studied the S = 1/2 spin chain system in which the interaction between 
neighboring spins changes as ferromagnetic- ferromag-netic-antifcrromagnetic (Fig.l). The Hamiltonian of 
this model is written by 

H = J A J2( h v-hv + Ah h-i,3j) 

3 

- J F X^i-Sj+l + Ah 3j,3j+1 
3 

+ ' l 3j+l,3j+2 + ^' l 3j+l,3j+2) i (1) 
fy,m = SfSm + > M,m = $1 &m > (2) 

where Ja > and Jf > are the magnitudes of the antifcrromagnctic and ferromagnetic couplings, respec- 
tively. I have introduced the anisotropy parameter A > for later convenience, although Hida investigated 
only the A = 1 case. This model is not too artificial nor a theoretical toy, as noticed by Hida. In fact, the 
substance 3CuCl2 • 2dioxane, of which magnetization is measured by Ajiro et al. 2 in strong magnetic fields, 
is known to be a (quasi-)one dimensional magnet to which this model is applied. 

Hida 1 investigated this model (only the A = 1 case) by the numerical diagonalization method for 
finite systems (up to 24 spins) to find that there was a plateau in the magnetization curve as far as 7 = Jf/ Ja 
is small, and the width of the plateau decreases as the parameter 7 increases (Fig. 2). The location of the 
plateau was Ms/3, where Ms is the saturation magnetization. He could not obtain a definite conclusion 
about the existence of the plateau for large 7 case, because the magnetization varied stepwise for the finite 
systems. Since this model becomes S = 3/2 antiferromagnetic chain model when A = 1 and 7 — > 00 due to 
the formation of an S = 3/2 quartet by spins S^j, S^j+i and S3J+2, it is expected that there is no plateau in 
this limit. It is believed that all the half-odd spin chains with antiferromagnetic nearest-neighbor interactions 
belong to the same universality class as that of S = 1/2. 

There arise several questions: 

(i) Why does the plateau appear at M/Ms = 1/3? (ii) Does the plateau vanish at a finite value of 7 or 
persist to 7 = 00? (iii) If the plateau vanishes at 7c, what is the behavior of the plateau near 7c? 

They are interesting questions not only from the standpoint of the statistical physics but also from that 
of explaining the properties of the existing materials. I note that the plateau was not observed in the report 
of Ajiro et al. 2 

In this paper I will analytically study the mechanism for the appearance of the plateau and also 
discuss whether the critical value 7c exists or not. 

Performing the spin rotation around the z-axis for the spins located at the sites 3j 

S%j —Sij ' &3j ^ ^^3] ' ^3j =^ ^fj j (3) 

Hamiltonian (1) is transformed into the form of a generalized version of the trimerized antiferromagne-tic 
chain: 

W=Jo(l-25i) J>£-_ li3j 
3 

+ Jb(l + S±) J2( h ksj+l + &«+l,3j+2) 

3 (4) 
+ J (A -25 z )Y^h z Sj _ h3j 

3 

+ Jo(A + 5 X ) J2( k h 3j+l + fc 3j+l,3j+2) . 
3 



2 



where 

2 J F + J A -2 7 +l A 

2 7 + 1 ' 2 7 + l 



In case of the XY model where A = 0, we can solve Hamiltonian (4) by the use of the fermion 
representation through the Jordan- Wigner transformation. The dispersion relation is shown in Fig. 3. When 
the magnetic field is applied, the dispersion curve shifts along the w-axis, which explains the existence of 
the plateau as well as its location M/M$ = 1/3. Note that the magnetization M is related to the number of 
occupied states as 

M 1 (occupied states) 



N 2 N 

where N is the total number of spins. 



(6) 



This is a simple explanation for the question why the plateau appears at M/Mg = 1/3. However, 
there still remains questions. Why is there no plateau in the experimental results on 3CuCl2 • 2dioxane? Is 
there any essential difference between the XY case (A = 0) and the isotropic Heisenberg case (A = 1)? How 
can we explain Hida's result, especially in 7 > 5 case? To answer these questions, we have to consider the 
interactions between fermions when A/0. 

To consider the effect of the trimerization and that of the interaction between fermions si- 
multaneously, I use the method of the bosonization. 3 The bosonization is one of the powerful methods in 
one-dimensional quantum problems. Here I do not enter into the details of the bosonization procedure. I 
note that the bosonization is usually done in case of no magnetization (half-filled in the language of the 
fermion), but in the present case, we have to perform the bosonization near M/M$ = 1/3 (1/3-filled in the 
language of the fermion) . 

After the bosonization, Hamiltonian (3) is transformed into a generalized sine-Gordon Hamiltonian: 

H = J dx{A(V6) 2 + CP 2 

-B ± cos6-B z (V6) 2 cos6} , (7) 

^).— ('-£)■ 

B ± = 26 ± /V3a 2 , B z = S z /tt , (9) 
[6{x), P(x')} =iS{x-x') , (10) 

where a is the distance between neighboring spins. The effect of the trimerization appears in the B± and 
B z terms. If the B± term and/or the B z term are relevant in the sense of the renormalization group, the 
spectrum of Tt has a gap, which brings about the plateau in the magnetization curve. If both of them are 
irrelevant, the spectrum of H is gapless, which results in no plateau. Which case is realized? — It depends 
on the magnitudes of B± and B z and also on the parameter 

77=27T- 1 V /CyI. (11) 

The renormalization group calculation shows that the B± term and/or the B z term are relevant as far as 
r\ < 4. 

The value of 77 is slightly shifted through the bosonization procedure. The expressions of A and B 
in eq.(8) is considered to be the lowest order expansions with respect to Aq. When Aq — » —1 (i.e. 7 — > 00), 
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the system becomes ferromagnetic. Therefore rj should diverge to +00 when Ao — > — 1, although it seems to 
diverge at A -> -tt/2\/3 = -0.907 from eq.(8). 

The bosonized Hamiltonian (7) has the same form as that of the generalized version of the 
dimcrizcd XX Z model. 



^J^{l + (-l)^i d) }^ +1 
1 

+ J^{A( d ) + (-l)^ d )}^ ;+1 , J>0, 
1 

In fact, if we perform the bosonization for Hamiltonian (12), we obtain 4 

fiW =Ja J dx{A^(\76) 2 + C^P 2 

- B^ cos 6 - B^ (V6») 2 cos 6} , 

where 



(12) 



(13) 



A^ = ±(l + ^) , cW=2*a(l-*?-) . (14) 

= S^/a 2 , B^ = . (15) 

Also in this case, the expression of and should be considered to be the lowest orde r expansio ns 
near A^ d ) = 0. In fact, in the absence of the dimerization, the exact form of r/ d ) = (2tt)^ 1 \/C^ /A( d ) is 
obtained from the application of the exact solution of the eight- vertex model as 5 



V 



(d) 



2/[l + (2/^)sin- 1 A (d) ] . (16) 



If we expand eq.(16) near A^ d ^ = 0, we can see that it agrees with the expression of r/ d ) from eq.(ll) and 
(14). 

TL has the same form as H^ d ' with the identification of the parameters 

A( d > = ^|, *i d) = ^|, W=6.. (18) 

Then we can use the knowledge on the dimerized XX Z model. The phase diagram 6-8 of the dimerized XX Z 
model when 5^ = 5^ = is shown in Fig. 4. This phase diagram was obtained by the renormalization 
group calculation, 6 by the high temperature series expansion after mapping onto the finite-temperature 
classical 2D model (modified Ashkin- Teller model), 7 and by the numerical diagonalization of the original spin 
Hamiltonian for finite systems. 8 When 5^ = 5i d \ a naive consideration leads to the effective dimerization 
parameter 

S^ = (28^+S^)/3, (18) 

because 5^ is related to S x and S y , whereas to S z . However, the renormalization group method and 
the variational method bring about 4 

5^ = (26^+r,6^)/(2 + r,) . (19) 

From eqs.(5), (17) and (19), we obtain the mapping of the present model onto the dimerized XXZ model, as 
shown in Fig. 4. The case 7 = 1 of the present model corresponds to the case A( d ) = -2/3V3 = -0.385, 5^ 
= 0.374 of the dimerized model, and the case 7 = 00 to the case A^ d ) = —1, 5^ = 1/2. Therefore there exists 
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the critical value 7c, where the transition from the plateau state to the no- plateau state. This transition is 
of the Kosterlitz-Thouless type, as known from the critical properties of the sine-Gordon Hamiltonian. 9 

The above discussion is based on the bosonization method, which make it difficult to estimate 
the value of 7c itself. It is because the parameters are slightly shifted through the bosonization procedure, 
as already explained. A rough estimation of the value of 7c is 

7c = 4 - 5 . (20) 

I have analytically investigated the plateau in the magnetization curve of the S = 1/2 fcrromagnctic- 
ferromagnetic-antiferromagnetic spin chain, which is first pointed out by Hida 1 by the use of the numerical 
diagonalization. The present analytical study semi-qualitatively explains the numerical result of Hida. 
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Figure Captions: 



Fig.l. Ferromagnetic-ferromagnctic-antifcrromagnetic spin chain. Solid lines and dotted lined represent 
ferromagnetic couplings Jf and antiferromagnetic couplings J a, respectively. Three spins in an ellipse make 
an S = 3/2 quartet when A = 1 and 7 — > 00. 

Fig. 2. Sketch of Hida's numerical result for the magnetization curve. The case (a) is for smaller 7 and 
(b) for larger 7. 

Fig. 3. Dispersion relation of H of cq.(4) in case of A = 0. (a): S > case, (b): S — case. 

Fig. 4. Mapping of the present model onto the dimerized XX Z chain. The open circle corresponds to 
the 7 = 1 case and the closed circle to the 7 = 00 case. 
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